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Abstract 

We study leptoproduction of J ftp at large Q 2 within the nonrelativistic 
QCD (NRQCD) factorization formalism. The cross section is dominated by 
color-octet terms that are of order a s . The color-singlet term, which is of or- 
der a?, is shown to be a small contribution to the total cross section. We also 
calculate the tree diagrams for color-octet production at order a 2 s in a region 
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of phase space where there is no leading color-octet contribution. We find 
that in this regime the color-singlet contribution dominates. We argue that 
non-perturbative corrections arising from diffractive leptoproduction, higher 
twist effects, and higher order terms in the NRQCD velocity expansion should 
be suppressed as Q 2 is increased. Therefore, the color-octet matrix elements 
(Og ( 1 5'o)) and (Of ( 3 -Po)) can be reliably extracted from this process. Fi- 
nally, we point out that an experimental measurement of the polarization of 
leptoproduced J/ip will provide an excellent test of the NRQCD factorization 
formalism. 
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I. INTRODUCTION 



Quarkonium production has been the focus of much experimental and theoretical at- 
tention in recent years. This surge in interest is largely due to the observation of gross 
discrepancies between experimental measurements of J/ if) and if)' production at the Collider 
Detector Facility (CDF) at the Fermilab Tevatron [JIJ and calculations based on the Color- 
Singlet Model ||. The dramatic failure of the traditional technique for calculating quarkonia 
production and decay rates has led to a new paradigm for understanding quarkonia: the 
nonrelativistic QCD (NRQCD) factorization formalism of Bodwin, Braaten, and Lepage |3J. 

A central result of the NRQCD factorization formalism is that inclusive quarkonium pro- 
duction cross sections have the form of a sum of products of short-distance coefficients and 
NRQCD matrix elements. The short-distance coefficients are associated with the production 
of a heavy quark-antiquark pair in specific color and angular-momentum states. They can 
be calculated using ordinary perturbative techniques. The NRQCD matrix elements param- 
eterize the effect of long-distance physics such as the hadronization of the quark-antiquark 
pair. These can be determined phenomenologically. 

The power of the NRQCD formalism stems from the fact that factorization formulas for 
observables are expansions in the small parameter v, where v is the average relative velocity 
of the heavy quark and anti-quark in the quarkonium bound state. For charmonium v 2 ~ 0.3, 
and for bottomonium v 2 ~ 0.1. NRQCD f-scaling rules Q allow us to estimate the relative 
sizes of various NRQCD matrix elements. This information, along with the dependence 
of the short- distance coefficients on a s and a, permits us to decide which terms must be 
retained in expressions for observables to reach a given level of accuracy. At low orders, 
factorization formulas involve only a few matrix elements, so several observables can be 
related by a small number of parameters. 

Prior to the innovations presented in Ref. 0, most J /if) production calculations took 
into account only the hadronization of cc pairs initially produced in a color-singlet 3 Si state, 
as parameterized by the NRQCD matrix element (Of ( 3 Si)). An important aspect of the 
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NRQCD formalism is that, in addition to the color-singlet contribution, it allows for the 
possibility that a cc pair produced in a color-octet state can evolve nonperturbatively into 
a J/ij} or ip'. This color-octet mechanism is central to the current theoretical understanding 
of charmonium production at the Tevatron. 

For a majority of phenomenological applications of the NRQCD factorization formalism, 
the most important color-octet matrix elements are (Of ( 3 «Si)), (Of ( 1 5 , )), and (Of ( 3 Pj)). 
They describe the non-perturbative evolution of a color-octet cc pair in either a 3 Si, 1 So, 
or 3 Pj angular momentum state into a J/if>. Using heavy quark spin symmetry relations 
0, it is possible to express all three P-wave matrix elements in terms of one: (Of ( 3 Pj)) = 
(2 J + 1)(C^( 3 P )) + 0(v 2 ). Thus, at leading order in v, there are three color-octet matrix 
elements. At this time they cannot be computed from first principles, and must therefore 
be extracted from experimental data. 

The most precise determination of the color-octet matrix element (Og ( 3 Si)) comes from 
a fit of a leading order theoretical calculation |5],|6|] to CDF data on J/ip production at high 
transverse momentum (P±). A linear combination of the remaining two color-octet ma- 
trix elements (Of^So)) and (O^( 3 P )) 

can also be determined by fitting the leading order 
calculation to CDF data on J/ip produced at low to moderate P±. The fit is extremely 
sensitive to theoretical uncertainties, so that the value determined for the linear combina- 
tion of (Cg ( 1 5'o)) and (Og ( 3 P )) can on ly be regarded as an order of magnitude estimate. 
There are other J/ip production processes, such as those measured at fixed target experi- 
ments in photonic JF], and hadronic || collisions, that could, in principle, provide a means 
to determine (Of ( 1 5'o)) and (Of ( 3 Po))- However, in each case there are large theoretical 
uncertainties that only allow for an order of magnitude estimate of these color-octet matrix 
elements. 

It is clear that testing the NRQCD factorization formalism requires a precise determi- 
nation of each of the leading matrix elements. Specifically, the challenge lies in an accurate 
measurement of (0$ ( 1 5'o)) and (0$ ( 3 Po))- in this paper we study leptoproduction of J/ip, 
and show that it is possible to determine these matrix elements from a measurement of the 
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unpolarized J/ip production cross section. We also show that once these matrix elements 
are extracted, the polarization of the J/ip can be predicted without introducing any new 
parameters. Thus, a measurement of the polarization of leptoproduced J / ip will provide an 
excellent test of the NRQCD factorization formalism. 

An explicit calculation need not be carried out to understand why it should be possible to 
determine (Of i}S )) and (Og ( 3 P )) from a measurement of the J/ijj leptoproduction cross 
section. The 0(a s ) contribution to J/ip production comes from processes in which the cc pair 
is produced at short-distances in either a color-octet 1 S'o or 3 Pj configuration. It is not until 
one includes the 0(a 2 ) contributions that it is possible for J/ip to be produced through the 
production at short distances of a cc pair in a color-singlet 3 S 1 state. The color-octet matrix 
elements (Of ( 1 5 , )) and (Of ( 3 -Po)) a re suppressed by v 3 and t> 4 , respectively, relative to the 
color-singlet matrix element. Since the perturbative coefficient for color-octet production 
is enhanced relative to the perturbative coefficient for color-singlet production by a factor 
of ~ vr/a s , we expect the color-octet contribution to be roughly v 3 n/a s ~ 2 times the 
color-singlet contribution. In fact, upon completing the calculation, we find that color-octet 
contribution is about 4 times the size of the leading order color-singlet contribution. 

Note that since Q 2 can be large, leptoproduction is a better process from which to extract 
{Og ( 1 <S'o)) an d (Of ( 3 Pq)) than photoproduction and low-energy hadroproduction. The latter 
two processes lack any large scale other than the charm quark mass, and consequently, 
perturbative corrections to leading order calculations are larger. In addition nonperturbative 
effects, such as higher twist corrections to the parton model and diffractive J / ip production, 
are less effectively suppressed in photoproduction and low-energy hadroproduction than in 
leptoproduction at large Q 2 . 

Finally, we point out that there are J/ip leptoproduction final states which can only be 
described by 0(a 2 ) tree-level contributions (i.e. there is no 0(a s ) color-octet contribution). 
The final states are those in which a gluon jet is well separated from the J/ip. Though the 
dominant contribution to these final states comes from the 0(a 2 ) color-singlet term, there 
are also 0(a 2 ) color-octet contributions. We calculate these color-octet terms to determine 
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if there is any region of phase space where they are enhanced relative to the color-singlet 
term. We find that in the regime where the calculation is valid the color-singlet term always 
dominates. 

This paper is organized in such a way as to separate theory from phenomenology. Readers 
who are only interested in the phenomenological implications of this work can skip the 
theoretical discussion without loss of coherence. 

We discuss theoretical issues in sections II-IV. First, we briefly review the NRQCD 
factorization formalism. Then we present the 0(a s ) and 0(a 2 ) tree level calculations of J/if> 
leptopro duct ion. At 0(a s ), we calculate both the polarized and unpolarized J/ijj production 
cross sections. Finally, we discuss possible corrections to our calculations from diffractive 
production and higher twist corrections to the parton model. We also address the possibility 
of a breakdown of the NRQCD velocity expansion near the boundaries of phase space. For 
the distributions computed in this paper, we argue that these errors are systematically 
reduced as Q 2 is made large. 

The phenomenology of J/ip leptoproduction is presented in sections V and VI. First, 
we discuss the present status of the determination of the leading order NRQCD matrix 
elements. Then we present our predictions, and discuss the theoretical limitations on the 
accuracy to which (OffSo)) and (Of{ 3 P )) can be measured in leptoproduction. We con- 
clude that leptoproduction of unpolarized J ftp is an excellent way to measure (G$ ( 1 5'o)) 
and (Of ( 3 P )), and that leptoproduction of polarized J/ip provides a powerful test of the 
NRQCD factorization formalism. 

II. THE NRQCD FACTORIZATION FORMALISM 

The NRQCD factorization formalism of Bodwin, Braaten, and Lepage || has emerged 
as a new paradigm for computing the production and decay rates of heavy quarkonia. This 
formalism provides a rigorous theoretical framework which systematically incorporates rel- 
ativists corrections and ensures the infrared safety of perturbative calculations |J. In the 
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NRQCD factorization formalism, cross sections for the production of a quarkonium state H 
are written as 

*(^ = £-^<°|G?l°>. (i) 

where tuq is the mass of the heavy quark Q. The short-distance coefficients, F n , are asso- 
ciated with the production, at distances of order l/rriQ or less, of a QQ pair with quantum 
numbers indexed by n (angular momentum 23+1 Lj and color 1 or 8). They are computable in 
perturbation theory. In Eq. ([]]), (0|O^|0) are vacuum matrix elements of NRQCD operators: 

(0|Of |0> = EE(°l4i^(A) +X)(H(X) +X\lC n \0}, (2) 
x x 

where K. n is a bilinear in heavy quark fields which creates a QQ pair in an angular-momentum 
and color configuration indexed by n. The bilinear combination K,\JC n has energy dimension 
d n . The production matrix elements describe the evolution of the QQ pair into a final state 
containing the quarkonium H plus additional hadrons (X) which are soft in the quarkonium 
rest frame. 

Throughout the remainder of the paper, we use a shorthand notation in which the vacuum 
matrix elements are written as (Ojl 8 ^( 2S+l Lj)). 

The NRQCD matrix elements obey simple scaling laws HQ with respect to v, the relative 
velocity of the Q and Q. Therefore, Eq. (|l|) is a double expansion in v and a s . Since the 
NRQCD matrix elements are sensitive only to large distance scales they are independent 
of the short- distance process in which the Q and Q are produced. Thus, it is possible to 
extract numerical values for the NRQCD matrix elements in one experiment and use them 
to predict production cross sections in other processes. 

III. LEPTOPRODUCTION CALCULATION 

At leading order in a s , J/ip are produced through the hadronization of a cc pair in either 
a 1 S'o, or 3 Pj configuration. The Feynman diagrams are shown in Fig. |TJ. 
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FIG. 1. Leading order diagrams for color-octet leptoproduction of J/ip 

Theoretical predictions for the J/ip leptoproduction cross section can be computed us- 
ing the techniques of Ref. [IIJ. The short- distance coefficients appearing in Eq. ([!]) can be 
computed by matching a perturbative calculation in full QCD with a corresponding per- 
turbative calculation in NRQCD. The production cross section of a cc pair with relative 
three-momentum q is computed in full QCD, and Taylor expanded in powers of q. In this 
Taylor expansion, the four-component Dirac spinors are expressed in terms of nonrelativis- 
tic two-component heavy quark spinors. The NRQCD matrix elements on the right-hand 
side of Eq. (|1|) are easily expressed in terms of the two-component heavy quark spinors and 
powers of q. The F n appearing in Eq. ([!]) are then chosen so that the full QCD calculations 
and NRQCD calculations agree. Detailed examples of the calculational technique can be 
found in Ref. ]TD] ; here we will merely quote our result for the differential cross section. 



The expression for the cross section determined from the diagrams in Fig. |I] is 



trie 



+ p^e + ij + X) = J J dy J 'dx f g/p (x) 5{xys - (2m c ) 2 - Q 2 ) 
2a s (^ 2 )a 2 e^ 2 1 1 + (1 - y f r . 3Q 2 + 7(2m c ) 2 ((^( 3 P )> 

]\ U 8 { b 0)) + 



m c (Q 2 + (2m c ) 2 ) { y L N Q 2 + (2m c ) 2 m 2 c 

_ 8(2m c ) 2 Q 2 (Oj( 3 P )) \ 
V (Q 2 + (2m c ) 2 ) 2 m 2 J ' U 

where s is the electron-proton center-of-mass energy squared, /i 2 = Q 2 + (2m c ) 2 is both the 
factorization and renormalization scale, and f g / p {x) is the gluon distribution function of the 
proton. The momentum fraction of the virtual photon relative to the incoming lepton is 
y = Pp-q/ P p -k, where P p is the proton four-momentum, q is the photon four-momentum, and 
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k is the incoming lepton four- momentum, and Q 2 = —q 2 . Note that the relative importance 
of (Of ( 1 5'o)) and (Of ( 3 -Po)) changes as a function of Q 2 . Thus it is possible to fit the 
differential cross section as a function of Q 2 and extract both of these matrix elements. 

The result presented in Eq. (|3]) holds for all values of Q 2 . Taking the limit Q 2 — > 
one recovers the photoproduction cross section |7| convoluted with the electron splitting 
function: 

a r dO 2 r 1 l + (l-v) 2 
lim a(e + P — e + V + X) -> — / \ dy —± ^a( 7 + P -> ^ + X) . (4) 



Q 2 ^o ' 2n J Q 2 Jo y 

As mentioned in the introduction, corrections to the photoproduction cross section from 
higher order perturbative QCD corrections terms, from diffractive photoproduction, and 
from higher twist effects, may all be large. However, in the high-energy limit Q 2 , s ^> (2m c ) 2 
we expect corrections to be negligible. Letting Q 2 ,s 3> (2m c ) 2 in Eq. (H) we obtain 

lim a(e + P ^ e + ip + X) -> [ I dy [ dx f g/N (x) 5(xys - Q 2 ) 

! 2 ,mi/s- j >0 J U J J 



m 2 /<3 2 ,m 2 /s-+0 J Q 2 



x 



2a s (Q 2 )a 2 e 2 c n 2 1 + (I - y) 2 ( {0 *,i s )) + 3 (g|gfi ))' 



Once the color-octet matrix elements have been extracted from measurements of the total 
cross section, they can then be used to make predictions for the polarization of J/ip produced 
in leptoproduction without introducing any new free parameters. The measurement of the 
polarization will be an important check of the NRQCD factorization formalism. The cross 
section for production of longitudinally polarized J/ip, where the polarization axis is the 
direction of the three-momentum of the J/ip in the photon-proton center-of-mass frame, is: 

dQ 2 



cr(e + p -> e + ip L + X) = J J dy J dx f g / p (x) 5(xys - (2m c ) 2 - Q 2 

\(OtCSo)) + 3- 



Q'' 

2a s ^ 2 )a 2 e 2 y j 1 + (1 - y ) 2 {/ ^ n g ^ , Q ((^( 3 P )> 
3m c (Q 2 + (2m c ) 2 ) { y 

, 1-y 48(2m c ) 2 Q 2 (Ot(*P )) 



m 2 



(6) 



y (Q 2 + (2m c ) 2 ) 2 m 2 J " 
The polarization can be measured by studying the angular distribution of the leptons in 
the leptonic decay of the J/ip. If 9 is defined to be the angle between the momentum of the 
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leptons in the J/ip rest frame and the momentum of the J /if) in the photon-proton center 
of momentum frame, then the decay distribution of the leptons is given by: 



dcos9 

where, 



oc 1 + a cos 6 (7) 



_ l-3/ L _ a(ip L ) ( , 

1 + /i cr(V') 

Note that the theoretical prediction depends only on one parameter, the ratio R = 
(e^( 3 Po))/(m2((9^(%))). In the limit where Q 2 ,s > (2m c ) 2 , the polarization parame- 
ter, a, goes to zero. This can be seen by inspection of Eqs. (|5|) and (|6|). 

At 0(a 2 ), J/i/j is produced through the hadronization of a cc pair in either a 3 Si, 1 S , 
or 3 Pj configuration . The Feynman diagrams are shown in Fig. |2|. Note that only the 
diagrams in Fig. |2]a produce a cc pair in a color-singlet 3 Si state; all other diagrams produce 
a color-octet cc pair. 

Once again we use the techniques of Ref. (H| to compute the J / ip production cross section 
from the Feynman diagrams shown in Fig The expressions obtained are complicated, so 
we do not present them hereQ However, we do report on checks we have made to ensure 
that that these expressions are correct. First, we have checked that in the Q 2 — > limit we 
recover the photoproduction cross section convoluted with the electron splitting function. 
Second, we numerically compared our color-singlet result to a calculation carried out by 



Merabet, Mathiot, and Mendez-Galain [O], and find agreement. 



IV. NON-PERTURBATIVE AND DIFFRACTIVE CONTRIBUTIONS 

Before presenting our predictions, we wish to discuss non-perturbative corrections to the 
NRQCD factorization formalism. Precise determination of the NRQCD matrix elements will 



1 The FORTRAN code for generating the differential cross sections presented later in this work is 
available by request. 
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(a) 



+ . . . 



(b) 




FIG. 2. 0(a 2 s ) diagrams for leptoproduction of J/ip. There are six diagrams of the type shown 
in (a), two diagrams of the type shown in (b), and two diagrams of the type shown in (c). The 
remaining quark diagrams, which are not shown, can be obtained from the diagrams in (b) by 
replacing the external gluon lines with quarks. Note, that only the diagrams in (a) contribute to 
the production of a cc pair in a color-singlet 3 Si configuration. 
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be impossible if these effects are not controlled. In the first part of this section, we consider 
diffractive leptoproduction and higher twist corrections to the parton model. In the second 
part, we discuss the breakdown of the NRQCD velocity expansion near the boundaries of 
phase space. 



A. Diffraction and higher twist 

At the HERA collider at DESY, diffractive processes contribute roughly 40% of the total 
7 — p cross section |12j| , and may be an equally important contribution to leptoproduction 
of J/ip. Moreover, the J/ip produced via this mechanism have similar kinematics to J/ip 
produced via the 0(a s ) color-octet contribution. 

To understand why diffractive leptoproduction of J / ip is kinematically similar to leading 
order color-octet production, it is necessary to introduce the following variable: 

Pv^± (g) 

Here is the four-momentum of the J/ip. In the proton rest frame z = E^/E^. For the 
process e + p — > e + J/ip + X, z = l + t/(s + Q 2 ), where we neglect the proton mass, and the 
mass of the final state system X. Diffractive processes are exponentially suppressed away 
from t = 0, where z = 1. 

The short- distance part of the 0(a s ) color-octet contribution leads to the production of a 
cc pair with no other particles in the final state; hence one would expect z = 1. However, soft 
gluons with momentum of order m c v 2 are radiated during the non-perturbative evolution of 
the cc into the J/i/j, resulting in z ~ 1 — v 2 . Thus there is significant overlap of the kinematic 
regimes of the 0(a s ) color-octet contribution and the diffractive contribution. It is possible 
to eliminate the diffractive contribution by requiring z C 1; however the 0(a s ) color-octet 
mechanism will also be eliminated by such a cut. There are other color-octet contributions 
at z <C 1, but they only occur at higher orders in a s where there is no longer a perturbative 
enhancement of the color-octet term relative to the color-singlet term. 
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Because diffractive leptopro duct ion is characterized by a large scale, Q 2 , this process 



can be studied using perturbative QCD. In Ref. ||13|| , a perturbative analysis of diffrac- 
tive leptoproduction predicts the cross section to fall as l/(Q 2 + (2m c ) 2 ) 3 , as compared to 
1/(Q 2 + (2m c ) 2 ) 2 for the color-octet mechanism. Therefore, at sufficiently large Q 2 , the 
diffractive contribution should be a negligible correction to our calculation. Furthermore, 



Ref. fl3| predicts that diffractively produced J / ip will be longitudinally polarized in the limit 



Q 2 » (2m c ) 2 . As discussed earlier in this paper, the 0(a s ) color-octet mechanism predicts 
leptoproduced J/ip to be unpolarized in this limit. Therefore, polarization of the J/ip will 
be a useful tool for distinguishing diffractive production from color-octet mechanisms. 

There are also obvious kinematic differences between color-octet and diffractive pro- 
duction which should make it easy to distinguish between the two. The most important 
distinction is that color-octet production will not lead to a rapidity gap in the final state, 
the hallmark signature of diffractive production. In color-octet production, the octet cc pair 
produced in the short-distance process must radiate at least one soft gluon before hadroniz- 
ing into the final state J/ip. This soft gluon, with energy in the quarkonium rest frame 
of order m c v 2 , will be well-separated in phase space from the proton remnant. A simple 
computation shows that the invariant mass of the soft gluon emitted by the color-octet cc 
combined with the proton remnant will be approximately vW 1P , where W 7P is the center- 
of-mass system energy of the virtual photon and proton. For the HI experiment at HERA, 
W 7P ranges between 30 GeV and 150 GeV ||14j| . If we take v 2 = 0.25, then for the color-octet 



contribution to e + p — > e + J /tp + X, we expect M x > 15 GeV. This is obviously much 
greater than Mx = 1 GeV, what is expected from elastic diffractive events. For diffractive 
disasso dative events, Mx can be greater than 1 GeV, but the cross section is expected to 
fall off as \jM\ |TJ]]. By requiring that there be no rapidity gap in the event and that 
M x 3> lGeV, it should be possible to extract a clean signal for color-octet J/ip production. 

In addition to the diffractive contribution, we must also consider the possible contribution 
from higher twist corrections to the parton model. These higher twist effects can be of great 
importance in quarkonium production calculations. 
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Well-known factorization theorems [J5j show that the inclusive production of a hadron 
can be written as a convolution of a parton distribution function, a parton scattering cross 
section calculable in perturbative QCD, and a fragmentation function which describes how 
the parton produced in the short distance process evolves into the final state hadron. This 
factorization formula receives higher twist corrections that scale as (KqcdI /•*)"> where /i 
is a characteristic scale of the process. For inclusive single hadron production, fi is P±, 
the momentum of the hadron transverse to the beam axis. By requiring that P± be large 
compared to Aqcd, we can ensure that predictions of the parton model calculation will be 
accurate. For quarkonium production, the situation is more complicated because there are 
two large scales that characterize the process: the mass of the quarkonium state, Mh ~ 
2rriQ, and Pj_. If some higher twist corrections are suppressed by Mu rather than Pj_, 
nonperturbative corrections to the parton model will not diminish as P± increases. This 
problem is particularly worrisome in the case of charmonium production, since 2m c rs 3 GeV, 
so nonperturbative corrections are not strongly suppressed. 

Higher twist corrections to color-singlet photoproduction are computed by Ma in Ref. flIS 



There it is shown that some of the twist-4 corrections are suppressed only by M|. Since 
higher twist structure functions are unknown, it is impossible to calculate the size of correc- 
tions, but it is clear that, at energies typical of the HERA collider, y/s 1P ~ 100 GeV, they 
can be important. Ma states in Ref. that in the case of electroproduction for fixed Q 2 
there are no twist-4 corrections that are suppressed by only Ml. However, no detailed cal- 
culation is given, and until such a calculation is carried out it is not clear what kinematical 
factor suppresses the higher twist contributions. We assume that higher twist corrections 
scale like Aq CD /(Q 2 + (2m c ) 2 ). An explicit calculation would be most welcome. 

B. Leptoproduction shape functions 

Another class of nonperturbative effects is associated with higher orders in the NRQCD 
velocity expansion |T7|. For quarkonium production near the boundaries of phase space, 
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it is sometimes the case that contributions from NRQCD operators which are higher order 
in v are enhanced by kinematical factors. This results in the breakdown of the NRQCD 
expansion. The crux of the problem is that in the perturbative QCD part of the matching 
calculation one uses twice the heavy quark mass instead of the quarkonium mass to compute 
the phase space for the production of the quarkonium meson. The difference between 2mq 
and Mh is a v 2 correction, which is ignored in leading-order calculations. However, at the 
boundaries of phase space this difference becomes important, and it is necessary to sum 
an infinite number of NRQCD matrix elements. This resummation leads to a universal 
distribution function called a shape function. Because of the universality of the shape 
functions, it may be possible in the future to test NRQCD by comparing shape functions 
extracted from different quarkonium production processes. However, in the present paper 
we are interested in extracting numerical values for the leading color-octet matrix elements. 
Therefore, we want to make sure that the effects of higher order terms in the NRQCD 
expansion are genuinely suppressed by powers of v 2 and can be safely neglected. 

The leading order contribution for producing J/ip through the hadronization of a color- 
octet cc pair in a 1 S configuration is: 



,A + (l-y)\( Q 2 + (2m c ) 



dxS(x, Q ) 5 y 



V \ xs j 

x J2(o\^T a x \J/^ + x)(j/ip + x\ x ^Ty\o) 

X 



where 



s , x Q2-, _ 8irf g/P (x)a s a 2 e 2 c 



9m c xsQ 2 (Q 2 + (2m c ) 2 ) ' 

The Ot( 3 P ) term can be analyzed in an analogous manner. In Eq. (fLOf) , = p c +p g , and 
in evaluating the phase space integral we have used the approximation = (2m c , 0, 0, 0) 
in the quarkonium rest frame. If we allow for the c and c produced in the short distance 
process to have nonvanishing velocity in the quarkonium rest frame, Eq. flip]) becomes: 
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da 



dQ 2 dy 



j y \ xs xs J 

x ^(0|^T a x | J/V> + X) (J/V» + X| X t T >|0) . 



In this expression, = (y^) + m c — m cJc(c)) is the charm (anti-charm) quark four- 
momentum in the J/ip rest frame, and A is a boost that takes us from the quarkonium 
rest frame to the frame in which we are performing the calculation. In order to perform the 
matching onto NRQCD, the delta function in Eq. (fLT|) is expanded in powers of 2P^-A(/ c +/ 5 ) 
and these factors are identified with the matrix elements of the NRQCD operators: 



xs 



n 



x 

Here, Dq = Do/(2m c ), where Do is the time component of the gauge covariant derivative. 

The infinite series of operators in Eq. can be summed into the universal shape function: 

do [, [, Q( ^ l + jl-yf j Q 2 + (2m c ) 2 8m 2 \ , 
dQ^ = J dX J dVESM } y 5 [ y xs 1^ VE ) F[ S ° 

(13) 



where 

F^S^Mve) = £(0|^ T a x |tf + X) (H + X\S (y E - iD ) X f T>|0). (14) 
x 

We now wish to consider the singly differential distribution da/dQ 2 . Integrating over y 
in Eq. flT3]), and expanding in powers of y E we find: 

d ° f, f, f, cv ^l + (l-2/) 2 z,n 0(8)1/ sj Q 2 + (2m c ) 2 8m 2 c ^ 



jL=JdxJdyJdy E S(x,Q 2 )—± y ± F ^S^](y E ) 5 ly 



■Ve 



dQ 2 J J J y \ xs xs 

= [dx S(x, Q 2 ) ( 1 + {1 ~ yLo)2 (O s CS )) - 2 -Zl^AF^ + — f](-)-A n F^) . 

J V VLO VLO VLO n=2 J 

(15) 

In this expression, y LO = (Q 2 + (2m c ) 2 )/xs, A = 2(2m c ) 2 / (Q 2 + (2m c ) 2 ), and 
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FW = Jdy E y E F?sP){y E )~v 2n . 



Notice that for Q 2 ^> (2m c ) 2 , A C 1 and the higher order corrections are suppressed by 
powers of Av 2 instead of v 2 . Thus the corrections associated with non-relativistic treatment 
of the quarkonium phase space can be made negligible by considering the singly differential 
distribution da/dQ 2 at large Q 2 . 

Higher order NRQCD corrections associated with the shape function can also be com- 
puted for the next-to-leading order color-octet and color-singlet contributions to leptopro- 
duction of J /if). The calculation is nearly identical to the calculation for photoproduction 



performed in Ref. fll? ], so we will merely quote results. The cross section is written in terms 
of the following variables: Q 2 , y, z, P±, and <f>. z is defined to be P p ■ P^/P p ■ q, where 
the non-relativistic approximation for P^ is used (P^ = (2m c , 0, 0, 0) in the J /if) restframe). 
This variable is distinct from experimentally observed z = P p ■ P^/P p ■ q, where P^ is the 
observed J /if} momentum. We define P± to be the momentum of the J /if) relative to the 
axis defined by the 7* three momentum in the 7*-proton center-of-momentum frame, and </> 
is the azimuthal angle of the J /if) about this axis. The differential cross section is given by: 



dQHydz J T J ± J 2H i zxs 

Lys - «•)(! - *) - _ Mi) 

V z [2m c )z[l - z) ) 

Eq. (|16|) applies equally well to color-octet and color-singlet production, so the quantum 
numbers of the cc produced in the short-distance process are not specified. It is worth 
noting that the shape function that appears in Eq. ( [L6| ) differs from the shape function in 
Eq. ([141): 5(yE — iD ) is replaced with 5(y+ — in - D), where is a light-like four vector. 
Expanding the delta function in Eq. (|T6| ) in powers of y + , it is easy to see that the moments 
of the shape function, = J dy + y r f_F[ 2S+1 L^j' 8S) ](y + ), are multiplied by a factor (1 — z)~ n . 
For 1 — z < v 2 , an infinite number of NRQCD matrix elements are relevant. For this 
reason, the differential distribution in z cannot be calculated in the forward region, even 
if Q 2 3> 4m 2 . However, other differential distributions, such as da/dQ 2 or da/dP]_, can 
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be computed because the integration over z provides sufficient smearing over the singular 
region. 

Based on the analysis presented in this section we have reason to believe that non- 
perturbative corrections are under control for leptoproduction of J/ip at large Q 2 . First, 
though the size of the diffractive contribution to leptoproduction is difficult to estimate 
because perturbative QCD calculations of this process are not expected to correctly yield 
the absolute normalization of the cross section [D|], perturbative QCD calculations do predict 
that diffractive production will be less significant at large Q 2 . In addition, rapidity gaps, 
the invariant mass squared of the final hadronic state, and the polarization of the J/ip 
produced in diffractive production, all differ in a qualitative way from color-octet production. 
These observables should allow experimentalists to isolate a clean color-octet signal. Next, 
higher twist effects are expected to scale as A.q CD /(Q 2 + (2m c ) 2 ). Finally, examination 
of higher order v 2 corrections associated with shape functions gives us insight into which 
differential distributions can be reliably computed within the NRQCD formalism. The 
distribution in z near the kinematic endpoint certainly cannot be predicted in the leading 
order approximation we consider in this paper. Therefore, we will not attempt to calculate 
da/dz, even though this quantity is often measured in experiments and has been considered 
by theorists in previous studies. The quantities which we do compute should be insensitive 
to these corrections. 



V. MATRIX ELEMENT PHENOMENOLOGY 



As discussed in the introduction, a naive counting of powers of v and a s tells us that 
the leading terms in the NRQCD factorization formula for the J/ip leptoproduction cross 
section are the 0(a s ) color-octet contribution, Eq. ([3]), and the 0(a 2 ) color-singlet contri- 
bution. The color-octet contribution depends on the NRQCD production matrix elements 
(Og^So)) and (Of ( 3 Pq)), and the color-singlet contribution depends on (Of ( 3 <Si)). These 
production matrix elements must be treated as parameters. Thus, before we present our 
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prediction for J / if) leptoproduction we discuss the current status of NRQCD J / if) production 
phenomenology. 

The color-singlet matrix element can be determined from the decay J /if) —>■ fi + n~ or 
from lattice calculations |18|]. The value determined for (Of( 3 Si)) is shown in Table |I[ 



The values of the color-octet matrix elements (Of ( 1 5 , )), (Of ( 3 P )), and (Og ( 3 Si)) have 
been extracted from CDF measurements [I] by Cho and Leibovich ||, and by Beneke and 
Kramer ||. Neither analysis extracts (Of ( 1 5'o)) and (Of ( 3 Po)) independently, instead they 
fit a linear combination, M k = {O^Sq)) + k (O 8 ( 3 ? ))/m c 2 , of the two. In Ref. f| k = 3, 
and in Ref. || k — 3.5. A detailed discussion of theoretical uncertainties in the extraction of 
these matrix elements is given in Ref. , and the results of their fit are shown in Table |. The 
largest uncertainty is due to scale setting ambiguities. Statistical errors are relatively small, 
only ~ 15% for (Os( 3 Si)), and ~ 25% for M 3 5 . In this analysis, the charm quark mass is 
chosen to be 1.5 GeV and errors associated with uncertainty in the choice of this parameter 
are not considered. Later we will consider errors in our calculation of leptoproduction of 
J /if), and we will find that uncertainty in the charm quark mass is a much greater source of 
error than uncertainty due to scale setting ambiguities. While this may not be the case in 
hadroproduction, we nevertheless feel that uncertainty in charm quark mass is a potentially 
large source of error which has not been included in the analysis of Ref. |§. This additional 
source of error must be kept in mind when considering the extraction of NRQCD matrix 
elements from this process. 

The extraction of M 3 5 appears to be particularly unreliable. The contributions to the 
cross section proportional to the matrix elements (Of ( 1 5'o)) an d (Of ( 3 -Po)) are subleading 
except at small transverse momentum, Pj_ < 5 GeV. Thus the extraction of these ma- 
trix elements is very sensitive to effects which can modify the shape of the P± spectrum. 
Specifically, the shape of the proton gluon density at small x can modify the slope of the 
P± distribution at low P±. This is reflected in the extreme sensitivity of the value of the 
parameter M 3 5 (see Table |) to the choice of the parton distribution function. The authors 
of take this as evidence that the result of the fit for M3.5 will be unstable to higher order 
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(Of ( 3 5!)) (GeV 3 ) 



(Ot( 3 S 1 )) (10- 2 GeV 3 ) 



M 3 . 5 (10-2 GeV 3 ) 



r(J/ip -> fi+fi-) 1.1±0.1 
CTEQ4L 
GRV(1994)LO 
MRS(R2) 

TABLE I. M3.5 is the linear combination (O 8 ( 1 S , )) + 3.5(Os( 3 Po))/m 2 c . The color-singlet 
matrix element is determined from the decay rate for J /if) — > The color-octet matrix 

elements are determined from a fit to CDF data performed in Ref. ||. First error is statistical, 
second is due to scale uncertainty. 

corrections and nonperturbative effects that will modify the P± distribution at low P±. As 
a result they conclude that (Og( l So)) and (Of ( 3 Po)) are not reliably extracted from CDF 
data. 

As discussed in the introduction (0% ( 1 5'o)) and (Of ( 3 Po)) can be extracted from ex- 
perimental data on photoproduction and low energy hadroproduction of J /if). However, 
theoretical uncertainties only allow for an order of magnitude estimate of the color-octet 
matrix elements. 



1.06±0.14t£gg 4.38±1.15I, 
1.12±0.14i^ 3.90±1.14l 
1.40±0.22lJ; 3 ^ 10.9±2.07± 



VI. RESULTS 

Our study of leptoproduction of J/ if) at large Q 2 will show that this process can provide 
a more accurate extraction of (Of ( 1 5'o)) and (Of ( 3 Po)) than hadroproduction or photo- 
production. Unfortunately, at this time, there does not exist appropriate data to which a 



reliable fit can be made. There does exist data on J /if) muoproduction [19|; a fit to the EMC 



data was carried out in Ref. ||20|| . However, we must caution that this data was taken with 



low to moderate values of Q 2 . In this regime, we expect large non-perturbative corrections 
to the NRQCD factorization formalism. More importantly, as we will see below, in the low 
Q 2 region the perturbative calculation suffers errors, due to the uncertainty in the charm 



20 



quark mass parameter. Therefore, we will not make a fit to the existing leptoproduction 
data. Instead, we will choose two sets of reasonable values for (0% ( 1 5 , )), and (Og( 3 P )) in 
order to make predictions. 

In the results presented below, we choose either (Of^So)) = 0.01 GeV 3 and 
(Og( 3 P ))/ml = 0.005 GeV 3 , inspired by the order-of-magnitude estimates presented in 
Table |, and in accordance with NRQCD f-scaling rules, or (Og ( 1 5 , )) = 0.04 GeV 3 , and 
{Ot( 3 P ))/m 2 c = -0.003 GeV 3 , which is consistent with the fit to low Q 2 EMC data car- 
ried out in Ref. ||20|| . In addition, we require Q 2 > 4 GeV 2 . This should reduce errors 
from higher order perturbative QCD terms and from higher twist effects. We also assume 
that experiments at HERA are most likely to make precision measurements of J/ip lep- 
toproduction, so we choose the center-of-mass energy to be y/s = 300 GeV, and require 
30 GeV < W lp < 150 GeV, where W 2 p = (P p + q) 2 . 

Given the above conditions, we calculate the 0(a s ) color-octet contribution to the J/ip 
leptoproduction cross section, and the 0(a 2 ) color-singlet contribution to the cross section. 
The results are presented in Table [TJ for two sets of kinematic cuts: Q 2 > 4 GeV 2 , and 
Q 2 > 4 GeV 2 , P\ > 4 GeV 2 , where P± is the transverse momentum of the J ftp in the 
HERA lab frame. The color-octet result is for the first (second) choice of (Og( l So)) and 
(Of ( 3 -Po)) given above. Note that with no P± cut the color-singlet contribution is ~ 1/4 the 
size of the color-octet contribution, while with a Px_ cut it is ~ 1/5 the size of the color-octet 
contribution. For either choice of color-octet matrix element values the 0(a 2 ) color-singlet 
contribution is small compared to the 0(a s ) color-octet contribution. 

If we consider final states in which a gluon jet is well separated from the J/ip, there is 
no contribution to the cross section from the 0(a s ) color-octet term. This is because the 
0(a s ) color-octet contribution produces J/ip through the emission of soft gluons from the 
cc pair produced in the hard scattering. These gluons have less than 1 GeV of momentum 
in the J/ip rest frame. However, there will still be a color-octet contribution from 0(a 2 ) 
tree level diagrams. To isolate events in which the gluon jet is well separated from the J/ift, 
we require that {P± P ) 2 > 2 GeV 2 and z < 0.8. P^ is the momentum of the J/ip transverse 
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Cuts 0(a s ) color-octet 0{a 2 s ) color-singlet 0{a 2 s ) color-octet 

89 pb 
40 pb 

13 pb 8 pb 

TABLE II. The 0(a s ) color-octet contribution to the J ftp leptoproduction cross section for 
(OfCSo)} = 0.01 GeV 2 and (^( 3 P ))K = 0.005 GeV 2 ((O$?S )) = 0.04 GeV 2 and, 
(Og ( 3 Po))/ml = 0.005 GeV 2 ), the 0{a 2 ) color-singlet contribution, and the O(o 2 ) color-octet 
contribution. The results are given for three different sets of kinematic cuts. Only the 0(a s ) 
color-octet calculation and the 0(a 2 ) color-singlet calculation are valid in the region of the first 
two sets of cuts. There are only 0(a 2 ) contributions in region of the third set of cuts. Note 
theoretical errors are not included. 

to the axis defined by photon three-momentum in the photon-proton center-of-momentum 
frame. The result of this calculation is shown in Table |I[ The 0(a 2 ) color-octet diagrams 
contribute about 40% of the total cross section. We do not show differential distributions in 
P]_, rapidity, and Q 2 , but report that the shape of these distributions is roughly the same 
for both 0(a 2 ) color-octet and 0{a 2 ) color-singlet contributions. Therefore, inclusion of 
color-octet mechanisms only serves to change the overall normalization of the total cross 
section in this region of phase space. Since the overall normalization is already uncertain in 
a leading order calculation, we do not regard events with J/ip and a gluon jet as particularly 
interesting for measuring color-octet matrix elements. 

In Fig. ||, we show the 0(a s ) color-octet contribution (upper histograms) and the 0(a 2 ) 
color-singlet contribution (lower histograms) to the differential cross section as a function 
of Q 2 . The solid histograms are without a P 2 cut, and the dashed histograms are with a 
P 2 cut of 4 GeV 2 . In generating the curves we used /i 2 = Q 2 + (2m c ) 2 , m c = 1.5 GeV, 
and (Oi^So)) = 0.01 GeV 3 , (Ot( 3 P )} /m 2 c = 0.005 GeV 3 . Next, we study the error in our 
calculation that results from uncertainty in the choice of renormalization and factorization 
scales. In Fig. |], we show the 0(a s ) color-octet differential cross section as a function of 
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Q 2 > 4 GeV 2 
Q 2 , Pi > 4 GeV 2 
(P7 P ) 2 > 2 GeV 2 , z < 0.8 



331 (226) pb 
290 (200) pb 



500 




10 20 30 40 50 

Q 2 (GeV 2 ) 



FIG. 3. 0(a s ) color-octet (upper histograms) and 0(a 2 ) color-singlet (lower histograms) con- 
tributions to the differential cross section as a function of Q 2 for two sets of kinematic cuts. The 
solid histograms were generated with a Q 2 cut of 4 GeV 2 , and the dashed histograms were generated 
cutting on both Q 2 > 4 GeV 2 , and P\ > 4 GeV 2 . 
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FIG. 4. 0(a s ) color-octet differential cross section as a function of Q 2 for three different choices 
of scale: fi 2 = Q 2 + (2m c ) 2 (solid histogram), /i 2 /4 (dotted histogram), and 4fi (dashed histogram). 

Q 2 for three choices of the renormalization and factorization scale: fi 2 = Q 2 + (2m c ) 2 (solid 
histogram), /i 2 /4 (dotted histogram), and 4/i 2 (dashed histogram). If we were to fit this 
result to experimental data the scale uncertainty would result in a theoretical uncertainty of 
about at Q 2 = 4 GeV 2 , and + 2 % at Q 2 = 50 GeV 2 in determining the values of the color- 
octet matrix elements. These small errors give us confidence that higher order perturbative 
corrections to the leading color-octet production mechanisms are under control. 

In Fig. |^, we study the error in our prediction for the 0(a s ) color-octet differential cross 
section resulting from uncertainty in the determination of the charm quark mass. The three 
histograms in the figure correspond to three different choices of m c : m c = 1.3 GeV (dotted 
histogram), m c =1.5 GeV (solid histogram), and m c = 1.7 GeV (dashed histogram). The 
error in our prediction for the cross section, as defined by 



is plotted in Fig. || for m c = 1.3 GeV (dashed line) and 1.7 GeV (dotted line). Though there 




24 




1 







10 



20 30 



40 



50 



(GeV" 



FIG. 5. 0(a s ) color-octet differential cross section as a function of Q 2 for three different choices 
of m c . m c = 1.3 GeV (dotted histogram), m c = 1.5 GeV (solid histogram), and m c = 1.7 GeV 
(dashed histogram). 

is some improvement as Q 2 is increased, the error is still ^5% a ^ Q 2 = ^ GeV 2 . This is 
much greater than errors we expect to receive from higher orders in perturbation theory and 
higher twist effects, so the uncertainty in the charm quark mass dominates the theoretical 
error in our prediction. 

In addition to the Q 2 distribution, we have also calculated the P]_ distribution shown in 
Fig. ^, and the rapidity distribution shown Fig. ||. The rapidity of the J/ip, is defined to be: 

y 2 \E-PJ' 

where E and P z are the energy and z-component of the J/ip three-momentum measured 
in the photon-proton center-of-momentum frame with the z-axis defined by the photon 
three-momentum in that frame. We only present results for one choice of scale, fi 2 = 
Q 2 + (2m c ) 2 , one choice of charm quark mass, m c = 1.5 GeV, and for (Of^So)) = 0.01 GeV 3 , 
{Oi^P ))/m 2 c = 0.005 GeV 3 . Note that in Fig. | there is no P x cut. The sharp fall off of 
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FIG. 6. Theoretical error due to uncertainty in charm quark mass as a function of Q 2 . Error 
is plotted for m c = 1.3 GeV and m c = 1.7 GeV. 

the distribution at P\ ~ 3 GeV 2 is due to the Q 2 > 4 GeV 2 cut. The P 2 and y distributions 
suffer considerable error due to the charm quark mass uncertainty. In the case of the Pj_ 
distribution, the errors are similar in size to the Q 2 distribution, and, like the Q 2 distribution, 
decrease slightly as P\ is increased. For the y distribution, the error is roughly 1^5% ■> ano ^ * s 
independent of y. 

Once the color-octet matrix elements have been determined from leptoproduction of 
unpolarized J/ip, it is possible to make definite predictions for the production of polarized 
J/i/j. The expression for the polarization parameter, a, is given in Eq. (Q) and depends only 
on the ratio of color-octet matrix elements R = (Of ( 3 Po)) /(m 2 (Og ( 1 >So))). In Fig. |9], we 
plot a as a function of Q 2 , choosing fi 2 = Q 2 + (2m c ) 2 , and m c = 1.5 GeV. Plots are shown 
for R = 1.0,0.5,0.25, —0.075. If the cc produced in the short distance process is in a x Sq 
state, the J/ip will be unpolarized, so we expect a = when R = 0. Positive values of R 
lead to slightly longitudinally polarized J ftp (a < 0), while if R < (i.e. (C^( 3 P )> < 0) the 
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FIG. 7. Differential cross section as a function of P±. Solid - 0(a s ) color-octet; dashed - 0{a 2 ) 
color-singlet. We require Q 2 > 4 GeV 2 ; no P± cut is imposed. 



J/ip emerges with slightly transverse polarization. The error due to the uncertainty in the 
charm quark mass is significant, but not as large as the errors in the P± or y distributions. 
To illustrate this, we plot the polarization in the case of R = 0.5 for three different values of 
m c in Fig. |10[ The error is largest at large values of Q 2 where it is approximately Note 
that we have only included the 0(a s ) color-octet contribution. Since the 0(a 2 ) color-singlet 
mechanism contributes about 25% of the total cross section, its effect on the polarization 
needs to be included before detailed comparison with experiment can be made. 

In summary, the total cross section for leptoproduction of J /if) is dominated by leading 
order color-octet terms; the 0(a 2 ) color-singlet term contributes only 20-25% of the total 
cross section. An accurate measurement of the distribution da/dQ 2 will, thus, allow for 
an independent determination of (Of ( 1 5'o)) and (Of ( 3 -Po))- As discussed in section IV, 
nonperturbative errors are expected to be small at large Q 2 . For Q 2 > 4 GeV 2 higher twist 
effects result in an error of less than 8%. In addition we estimate the error due to the scale 
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FIG. 8. Differential cross section as a function of rapidity (y). Solid - 0(a s ) color-octet; dashed 
- 0(a 2 s ) color-singlet. We impose the cut Q 2 > 4 GeV 2 

uncertainty; we find it is less than 10%. Finally we estimate the error due to the uncertainty 
in the determination of the charm quark mass, and find it to be the dominant error. One 
must require Q 2 > 10 GeV 2 for the error to be less than 50%. Errors in the extraction 
of (Og ( 1 5'o)) and (Og ( 3 -Po)) from photoproduction (Q 2 = 0) are obviously greater than 
100%. The errors that arise in extractions from hadropro duct ion reported in Ref. || are 
greater than 100%, and do not even include error due to uncertainty in the charm quark 
mass. Therefore, measurement of leptoproduction of J/ip at Q 2 > 10 GeV 2 will be a great 
improvement over existing extractions. Nevertheless, it is somewhat disappointing that our 
prediction is dominated by the error due to the uncertainty in the charm quark mass since 
it cannot be systematically reduced by doing higher order perturbative QCD calculations. 

Once (Cg ( 1 5'o)) and (Og ( 3 P )) are determined, distributions in P 2 and rapidity, as well 
as the polarization of the J /if}, can also be predicted. The polarization of J /if} is least 
sensitive to errors associated with the uncertainty in the charm quark mass. Thus, an 
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FIG. 9. Polarization parameter a as function of Q 2 . Dashed line - R = 0.25; Solid line - R = 
0.5; Dotted line - R = 1.0; Dashed-dotted line - R = -0.075 

experimental measurement of the polarization will be an important test of the NRQCD 
factorization formalism. 

In addition, we have shown that the production of J/ip with a well separated gluon jet 
occurs only at 0(a 2 s ), and makes a small contribution to the total cross section. For these 
final states, the color-octet contribution is suppressed, and only modifies the normalization, 
but not the shapes of distributions. For this reason, we do not feel that production of J/ip 
in association with a hard jet will shed any light on color-octet mechanisms. 

VII. CONCLUSION 

In this paper, we studied leptoproduction of J/ip at large Q 2 , and found that the cross 
section is dominated by 0(a s ) color-octet production mechanisms. Thus, by measuring the 
total cross section as a function of Q 2 , it is possible to obtain measurements of the color- 
octet matrix elements (Of^So)) and (Of ( 3 -Po))- Neither of these has been determined to 
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FIG. 10. Polarization parameter a as function of Q 2 . R = 0.5; Solid line - m c = 1.5GeV; 
Dotted line - m c = 1.7GeV; Dashed line - m c = 1.3GeV. 

better than an order of magnitude. The dominant error is due to the uncertainty in the 
value of the charm quark mass. This results in roughly a 50% uncertainty in the prediction 
for the differential cross section for Q 2 > 10 GeV 2 . We, therefore, estimate that (OfO-Sb)) 
and (Of( 3 P )) can be measured in J/vf) leptoproduction at about the 50% level. 

In addition, we pointed out that once the color-octet matrix elements have been measured 
it is possible to make a parameter free prediction of the polarization of leptoproduced J/ip. 
This will provide an important test of the NRQCD factorization formalism. 

As a result of the work presented in this paper we conclude that leptoproduction of un- 
polarized J/ip at Q 2 > 10 GeV 2 will provide a much better determination of the color-octet 
matrix elements {Of ( l S )} and (O$( 3 P )) than currently available. In addition a measure- 
ment of the polarization of leptoproduced J/ip will provide an important test of the NRQCD 
factorization formalism. We would, therefore, encourage an experimental measurement of 
J/ip leptoproduction at large Q 2 . 
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